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Abstract 

The ground state energy shifts and excitation spectra of charged 
buckminsterfullerene Cqq , n = 1, ... 5 are calculated. The electron- 
vibron Hamiltonian of Part I is extended to include all A g and H g 
modes with experimentally determined frequencies and theoretically 
estimated coupling constants. Complex splitting patterns of H g vibra- 
tional levels are found. Our results are relevant to EPR measurements 
of spin splittings in CgQ and CgQ in solution. Spectroscopic gas-phase 
experiments will be of interest for further testing of this theory. As 
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found in Part I, degeneracies in the electron and vibron Hamiltoni- 
ans give rise to a dynamical Jahn Teller effect, and to a considerable 
enhancement of the electronic pairing interaction. This helps to over- 
come repulsive Coulomb interactions and has important implications 
for superconductivity in K3C60 and the insulating state in K^Cqq. 

PACS: 31.30,33.10.Lb,33.20.E,33.20.F,71.38.+i,74.20.-z,74.70.W 

1 Introduction 

This paper continues the investigations of the electron-vibron interactions 
within a fullerene molecular anion, Cqq . In Part I (Ref. [Q) we considered 
in great detail the idealized case of a single H g vibron mode coupled to 
the electronic degenerate t± u orbital for n=l,..5 electrons. By solving the 
problem for strong, intermediate and weak coupling regimes, we have shown 
the existence of Berry phases for odd n, and their importance in determining 
ground state energies and degeneracies. The Berry phase effects are clearest 
in strong coupling, where the parity of the pseudorotational orbital angular 
momentum L, is related to the electron filling n by 



The effects on the energies however are relatively stronger in weak coupling 
where quantum corrections enhance pair binding which is a factor 5/2 larger 
than the corresponding classical Jahn- Teller (JT) relaxation energy E JT . 
This leads to larger electron-vibron pairing interaction than previously cal- 
culated using Migdal-Eliashberg theory. Moreover, this enhancement is of 
direct importance to C$q, where the electron-vibron coupling is weak to in- 
termediate. These encouraging results clearly call for a more realistic study 
of the full electron-vibron problem of CgJ". This is the purpose of this paper, 
where we will address both the vibronic spectrum, and the pairing energies, 
in quantitative detail. 




(1) 
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The full molecular Hamiltonian describes the dynamics of 60 carbon 
atoms plus 240+n valence electrons. The problem is substantially simpli- 
fied by assuming the knowledge of the non interacting spectra for electronic 
levels H H], £|] and for molecular vibrations, both from theory || || and exper- 
iment H ||. For n—0, the electrons form a closed shell. This allows a Born- 
Oppenheimer decoupling of the vibron and electron systems. For n=l,...5, 
however, the extra electrons partly fill the threefold degenerate Lowest Un- 
occupied Molecular Orbitals (LUMO) states of t\ u symmetry, which gives 
rise to a linearly coupled JT system. Other orbitals, such as the Highest 
Occupied Molecular Orbitals (HOMO) states of h u symmetry, at ~ — 2eV 
below the LUMO states, and the t2 9 (LUMO+1) states at ~ +leV, intro- 
duce weaker quadratic electron-vibron couplings which we presently ignore. 
Further complications may arise from anharmonic effects. 

Fortunately, neglect of all these higher order effects is expected to be a 
very good approximation in Cqq, where the t\ u level is very well separated 
from others, and the C-C bonds are rather stiff and harmonic. Detailed 
Hartree-Fock calculations have shown, for example, that the energy gain by 
going from I h to D 5d , D 3d , D 2 h symmetries via static JT distortions, are 
in fact identical to within 1% [0. Therefore, restricting to the t\ u orbital, 
to linear JT coupling and harmonic vibrons is very well justified in C$q. 
Neglecting also Coulomb interactions (they will be discussed separately in 
Section |3|), the full electronic problem is therefore replaced by a 3x3 matrix 
linearly coupled to vibrons. 

Thus our Hamiltonian is an extension of the single mode mode solved 
in full in Part I. Here we shall include the eight H g vibron modes of real 
Cqq, instead of only one. Also we shall include two A g vibrons which also 
couple linearly to the LUMO electrons, even though they do not split its 
degeneracy. Symmetry prevents all other vibrons different from A g and H g 
to interact linearly with the t\ u orbitals. If we further neglect higher order 
interactions, all other vibrons are decoupled and unaffected by changing the 
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electronic filling n. 

Generally speaking, even with these drastic approximations, a realistic 
description of the dynamical JT state of a Cqq ion is a more complicated af- 
fair than the single mode treated in Part I. For a general coupling strength, 
there is in fact no linear superposition between effects produced by different 
H g vibrons coupled to the same t\ u orbital. Luckily, however, linear superpo- 
sition turns out to be valid in the weak coupling perturbative regime, which, 
in turn, applies, even if only approximately, to fullerene. 

In this paper we apply perturbation theory to the full electron-vibron 
problem of CgJ", with n—1 to 5 in a t\ u orbital, and including all A g and 
H g modes ||. As in Part I, we make the further approximation of replacing 
a spherically symmetrical coupling to the true icosahedral hamiltonian on 
the Ceo ball. Analytical expressions correct to second order in the electron- 
vibron coupling strengths are found for ground state and excitation energies, 
as well as for electron pair binding energies ("effective Hubbard t/'s"). Using 
realistic coupling constants from Local Density Approximation calculations, 
we obtain numerical estimates for these vibron-induced pair energies and 
find them to be unexpectedly large and negative (~ — 0.2eV) for n=l,3,5, 
and even larger but positive (~ 0.4eV) for n=2,4. This finding is discussed 
in qualitative connection with superconductivity in A3C60 H HH]) an d 
with the insulating state of A4C60 (A=K,Rb). 

The calculated vibron spectrum of Cgg" is also presented in detail. Not- 
withstanding the uncertainty in the physical coupling constants, large split- 
tings of all H g modes are predicted. These are expected to be observable for 
example in gas phase Cg~ and ions. This part is organized as follows: 
Section | defines the multi vibron model. Section ^| describes the perturba- 
tive calculation of the spectrum. Section |4] presents the relation of the results 
to experimental measurements of vibron spectroscopy of Cg^ anions. Section 
|5| discusses the interplay between electron-vibron and Coulomb interactions 
for the ground state and pair binding energies. We conclude with a short 
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summary. 



2 The Hamiltonian 

The single electron LUMO states of Cgo are in a triplet of t\ u representation. 
The important vibrational modes which couple to this electronic shell are of 
two representations: A g (one dimensional) and H g (five dimensional). A g , t\ u 
and H g are the icosahedral group counterparts of the spherical harmonics 1, 
{Xim\ln=-ii an d {^2m}m=-2 respectively. By replacing the truncated icosa- 
hedron (soccer ball) by a sphere, we ignore lattice corrugation effects which 
are expected to be small for the electron-vibron interactions, since they do 
not lift the degeneracies of L = 0, 1, 2 representations. 
The Hamiltonian includes the terms 

H =H° +H e ~ v + •■• (2) 

where electron-electron interactions, anharmonic interactions between phonon 
modes, and anharmonic coupling terms have been neglected. The non inter- 
acting Hamiltonian is 

H° = hJ2^k ( b kM b kM + l) +(e- n)J2 c L c ms, (3) 

k M=-L k ,L k V A ' ms 

b kM creates a vibron of mode k and energy Uk in the spherical harmonic 
state Yh k Mi where denotes the angular momentum of mode k, either or 
2 according to whether k is an A g or an H g representation respectively, 
creates an electron of spin s in an orbital Y\ m - This Hamiltonian operates 
on the basis 

II \ n kAl)v II \n mS )e (4) 
KM ms 

where \nkM)v (\n m s)e) is a vibron (electron) Fock state. By setting /i — > e we 
discard the second term in (H). 
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The electron-vibron interaction is local and we assume it to be rotation- 
ally invariant. The nuclear vibration field of eigenvector k is 

= E -^(YL kM m b L k Mk + Y LkM {(l)b Lk Mk) (5) 

M V Z 

where Cl is a unit vector on the sphere. The interaction between the vibrations 
and electron density is 

H e ~» oc$> fc fdn Uk (n)J2^)M^) (6) 

k s 

where the electron field operators are 

1 



V> a (fl) = Y lm{tyCms (7) 

m=— 1 

Using the relation 

dQ Y LM (Q)Y lmi (Q)Y lm2 (Q) oc (-1) M (L,-M\l mi ;lm 2 ), (8) 



where (• • •) is a Clebsch-Gordan coefficient (T^j, yields the second quantized 
Hamiltonian 



H e ~ v = H e A - v + H e H - 



H x v = \ ^T,9^ k Y,(-^r(bi +b k0 

» Z fc=l ms 

x (0, 0|1, — m; 1, m)c] ns c ri 



V3 1,1 
2 



#2? = V*E^ E E(-!) m Km + -m 



fc=3 Af=-2 ^ 

x (2,M\l,-m;l,m + M)cl ns c m+M s , (9) 

where the numerical constants are fixed by the requirement that the classical 
JT energy gain of a single mode k is gl/2 
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3 The Multi-Mode Spectrum (Weak Coupling) 

The perturbation hamiltonian (HD written in the Fock basis (§J), connects 
states whose number of vibrons N v (k) of mode k differs by exactly ±1. Thus, 
first order corrections to the energies vanish. Second order corrections are 



obtained by diagonalizing the matrix [13 



Ag = (a\H^ E ^ l _ H H^\b) : (10) 

where \a) and \b) are members of the same degenerate manifold, i.e. they have 
the same number of vibrons N v (k). The sum implied by the inverse operator 
(E^ — if ) -1 extends just to those states whose N^s differ only by ±1 from 
that of the multiplet being perturbed. This means that to second order in 
the coupling constants there are no direct inter-mode interactions, and the 
modes can be treated separately. The only second-order inter-mode coupling 
is a consequence of all the modes having a common ground state (N v =0). The 
effect of vibron k affects this N v =0 state either with a pure shift or through 
both shift and splitting. However all other modes k', having their ladder 
built on the same N v =0 state, are shifted or split by vibron k, according to 
the same structure of this N v =0 multiplet. This effect takes place through 
additive contributions proportional to g\ to the diagonal matrix elements A^j, 
relative to the k! ladder, without involving off-diagonal inter-mode couplings, 
which would be related to g 4 and higher order corrections. 

A single A g mode coupled to a t\ u level is the simple polaron problem, 
which is exactly soluble fL4| : the second order energy is exact. Since the A g 



representation is one-dimensional, it does not split the electronic degeneracy. 
The only effect is a downward shift of the whole spectrum. The amount of 
Ag-related energy shift is found to be —Ej[ T , —AE( T and —6E£ T , for n—1, 
2 and 3 respectively (where E( T is the classical JT energy gain = g^huj^jl 
of that mode A g {k)). Obviously, these results hold for both unpolarized and 
polarized spin states, since the t lu levels remain degenerate. 
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For the H g modes the situation is more complicated. For n—1, the degen- 
erate vibronic t\ u ground state is not split by JT coupling, as it conserves its 
L=l symmetry. For this reason, the only contribution of the vibron H g (k) 
to the spectrum of another vibron H g (k') is just a constant energy shift of 
-5/4 g\, which obviously does not affect energy differences in the spectrum 
of vibron k'. This is perfectly analogous to the effect of a single A g mode on 
all the other vibrons. 

By contrast, the n=2,4 and n=3 lowest vibronic multiplets (N v =0) are 
split (into 1 S, X D, 3 P, and 2 P, 2 D, 4 S levels respectively). Correspondingly, 
the levels of an interacting vibron k receive different diagonal contributions of 
order g\, from different interacting vibrons k', giving rise to a more intricate 
pattern of splittings. 

For the sake of simplicity, and also since Cg seems easiest to obtain in the 
gas phase fi~5f , we will concentrate on the many-modes spectrum for n—1. 
The perturbative results for all H g and A g modes will be presented in Section 

In principle, the full spectra for n— 2,3,4, can be determined following the 
same method. 

The spectrum for a single H g mode is given in Table I for the degenerate 
multiplets N v =0 and N v —1, and for n—1, 2, (5=0) and 3 (5=1/2) electrons. 
In Figures ([l], ^ |||) we replot on an expanded scale the results of exact 
diagonalization of Ref [p]] for 1, 2 and 3 electrons, along with the straight 
lines corresponding to present perturbative results for the lowest few states. 
These figures confirm that perturbative results retain quantitative validity 
up to g ~ 0.3. In the special case n—1, moreover, the perturbative results 
lay within 0.05huj of the exact value up to g < 0.4. 

We have therefore an approximate analytical estimate of the splittings 
induced by JT coupling, valid in weak coupling. For example, for n—1 the 
H g vibron excitation, originally at energy 1 above the ground state, splits 
into three vibronic levels with relative shifts — fg 2 , §g 2 and |g 2 (in units 
of hu of that vibron). Table I contains the complete list of these low lying 
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excitations energies, accurate to order g 2 . 

Because the effects of all A g and H g modes can be linearly superposed, 
there are two ingredients only, which we need to have in order to transform 
the analytical shifts of Table I into actual numbers for Ceo~ : the frequencies 
huJk and the coupling constants g^ of each individual H g and A g mode. For 
the frequencies, there are both calculated and measured values. We can 
avoid uncertainties by choosing the latter as given, e.g. for neutral Ceo i n 
Ref. ||. In doing so, we neglect the well known small systematic frequency 
shifts associated with bond-length readjustments and other electronic effects 
going from Ceo to Ceo" HH, |i~7j . They also depend on the environment of the 



Ceo ion. 

There are several calculated sets of coupling constants g^ || [10], |9], [18| , but 



no direct measurement. Since the agreement among the different calculations 
is far from good, in Section f| we present the results of several selected sets, 
which provides an estimate of the relevant uncertainties. We will eventually 
adopt the most recent values of Antropov et al. ]TB| . 

As seen in Table II, almost all of the couplings are weak, g^ < 0.4. As 
discussed, in this range the perturbative results are accurate within ten per 
cent or better for all the low lying states. As the discrepancies among the 
various estimates of the gi is much larger, these perturbative formulas are at 
this stage more than adequate, and particularly good in the n=l case, where 
the distortion is smallest. 

This is fortunate, since exact diagonalization is computationally rather 
demanding if all H g modes are included. Better knowledge of frequencies 
and coupling constants might warrant a calculation of higher orders in g^. 

4 Vibron spectroscopy of Cg ~ anions. 

The electron affinity of Ceo is large (2.7 eV) and experimental evidence has 
been found that the Ce~ (T^j and C^o 13 are stable ions in vacuum. In 
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solution a wider spectrum of ionization states has been demonstrated elec- 
trochemically, up to and including C|o~ [0, 21, j22|. As an adsorbate on a 



metal surface, the electronegative Cgo molecule naturally picks up electrons 



23|, |24| , and recent evidence has been provided of charge transfer which can 



be as large as n=6 fl7fl . In the solid state, finally, there are compounds, 
covering a wide range of charge transfers, from n=l, as in TDAE + -Cg |25| 
or RbiCgo [HI, n—3, as in K 3 Ceo or RbsCgo |p7j , n=4 as in K 4 Ceo [p8| , n=6 
as in Rb 6 C 60 
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33 in. 



or even higher as in Li 12 Cgo 
Among these systems, our calculations so far address concretely only the 
gas phase case. Unfortunately, to our knowledge no investigation appears 
to have been made of the vibrational excitations of isolated Cg and C|q . 
Our calculated excitation spectrum for Cg therefore constitutes a prediction 
which we hope will stimulate new work. 

In Table II we report the excitation energies predicted by perturbation 
theory, applied to the eight modes in the Cg case. Selection rules are not 
discussed here for any particular spectroscopy. We simply give the symmetry 
assignments. 

As is seen, the predicted splittings due to dynamical JT coupling are 
generally quite large, and should be well observed spectroscopically. However, 
as indicated by comparison between different sets of g^'s, there is a large 
uncertainty in these predicted splittings of the same order of magnitude as 
the splittings themselves. As remarked earlier, the same uncertainty does not 
affect the energetics of the following section, which is on safer grounds. Our 
calculated spectrum is therefore of qualitative value, and we rather expect it 
to work backwards. That is, a future precise measurement of the splittings 
should provide an accurate evaluation of the actual couplings. 

As a further caution, we should stress that our spherical representations 
in the Hamiltonian (|9]) neglect interactions due to the icosahedral lattice of 
carbons. For example the vibron multiplet of L = 3 decomposes due to the 



lattice into T- 



etc. In addition to neglecting lattice effects and 
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anharmonic interactions, we also ignore spin-orbit coupling. As remarked in 
Part I, it has been shown j3!| to yield splittings of the order of 50 cm -1 to the 
L =1 ground state (in Ar matrix), which is not a negligible amount. Thus we 
estimate that the splittings obtained by our Hamiltonian should dominate 
the splittings found in the real spectrum. 



As pointed out by Bergomi and Jolicoeur ||33|| , experiments on anions in 
matrix may be relevant to the vibronic effects. Near infrared and optical 
spectra of Cg^" ions in solution are available | 21| . A major t lu — > t\ g optical 
transition near 1 eV is present for all n values. It is accompanied by additional 
vibronic shake-up structures, typically near 350, 750, 1400 and 1600 cm -1 . 
This limited information seems as yet insufficient for any relevant comparison 
with our calculations. Well defined vibrational spectra are instead available 
for chemisorbed Cg^ [0 and for A n Cgo alkali fullerides p9| . In this case, 
however, interaction of the electronic t\ u level with surface states or with 
other t\ u states of neighbouring balls must turn the level into a broad band, 
and our treatment as it stands is invalid. One can generally expect rapid 
electron hopping from a molecule to another to interfere substantially with 
the dynamical JT process, in a way which is not known at present. The 
spectra of charged C 60 adsorbates and solids, in any case, do not present 
evidence of any dynamic splittings such as those of Table II, but rather of 
a gradual continuous shift most likely due to a gradual overall change of 
geometry, as suggested also by LDA calculations ||. 

Summarizing, we are yet unaware of detailed spectroscopic confirmation 
of the electron-vibron effects. We expect however these effects to be observ- 
able in the gas phase of Cg and C|o • In particular for Cg~ an observation 
that the lowest H g vibron splits into a near 7-fold degenerate multiplet (i.e. 
from the L = 3 = T2 U © G u pseudorotation level) would be an important 
confirmation of the electron vibron theory and effects of Berry phases. 
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5 Ground state energetics and effective Hub- 
bard C/'s. 



Perturbation theory allows us to write analytic expressions for the energy 
gain of the ground state at different n, and therefore for the pair energy 

U n = E n+ i + E n -i — 2E n , (11) 

as discussed in Ref JJj. Comparison with exact single- mode results shows 
a systematical perturbative overestimate (Figures [1], ||, [3D of the if 5 -related 
ground state energy shift. The error is however relatively small and quite 
acceptable for the couplings in Table II. The shift is 5/2 times larger than 
its classical value, as discussed in Ref. [|I|]. This factor 5/2 is important, 
because it leads in turn to a surprisingly large energetic lowering even for 
small g's, making JT vibronic coupling a much more important affair than 
it was understood so far. The physical reason for the large energy gain is 
that the dynamically JT distorted molecule undergoes a dramatic decrease 
of vibrational zero-point energy. This adds an extra —^g 2 hu (in the n=l 
case, say) to the static JT gain —^g 2 hu of each H g mode. The zero-point 
energy decreases faster at small g probably because the mexican-hat potential 
well is more "square- well" -like than the original harmonic potential. We also 
note that the proliferation of excited states upon coupling H g with t\ u is of 
fermionic origin, and does not add to the zero-point energy. 

Within second order perturbation theory the ground state energy is a 
sum of all the 2+8 contributions of the A g + H g modes: 

2 10 

E tot (n) = E Ag + E Hg = a n E J A T g + b n E J H T g = a n £ E J k T + b n J^ E J k T , (12) 

k=l k=3 

where 

E( T EE Iflgfc*, (13) 
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We already discussed in Section |3| the JT ground state energy gains due to an 
A g mode, with coefficients a n —-l, -4, and -6, for n=l, 2 and 3 respectively. 
Table I gives the corresponding energies for the H g {k) modes (k > 2). These 
are given by 

-\Ef k n=l 



n 

15 



lOE^ k n = 2 



- Y E 2k n = 3 (14) 

the appropriate coefficients being therefore 6i=-5/2, &2=-10, &3=-15/2. 

These expressions allow us to compute the individual contribution to the 
pair energies U n (Eq. ([□])) due to the A g and H g modes. We give these 
formulae in Table III. The corresponding numerical values are reported in 
Table IV, based on the ground state energy gains as given by the set of 
coupling constants of Eq. ([TT|) . 

We consider the unpolarized spin sector. Similar subtractions could easily 
be done, if needed, for high-spin states, or high- and low-spin, using, for 
example E tot (n=2,S=0), with E tot (n=3,S=3/2). Although, as we pointed 
out, the values of the individual g^s of Ref. |L8), [|l(| and are significantly 
uncertain, the global E Hg is much less author-dependent, amounting to 102 
meV, 84 meV and 78 meV respectively (n=l). 

As Table II shows, the coupling with the A g mode pushes JJ\ further 
towards negative values, but has the opposite effect on U2 and U3. 

The overall enhancement factor 5/2 in the ground state H g shift ends up 
producing a much larger pair energy than expected so far based on classical 
JT energies 0. In particular, our calculated JT energy gain of ~ 0.4 eV for 
n=2 and w 0.3 eV for n=3 (low-spin) is almost one order of magnitude larger 
than the currently accepted values! This has important implications, first of 
all, in determining whether the simple CgJ" ion, in vacuum, in a matrix or in 
solution, will choose to be high-spin or low-spin. 
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In order to discuss this point, we recall the existence of an intra-ball 
Coulomb repulsion U (not to be confused with the pair energy of Ref 
which for a t± u level is a matrix, specified by two main values, Uu (two elec- 
trons in the same orbital), and U± < Un (two electrons in different orbitals). 
Using the JT energy differences between the high- and low-spin states as 

4Ef - -E J H T ) 

Ag 2 H 3 ' 

u [ l ] ) -0.3ey (15) 

- (3C/f - 6££) 

t/| 3) ) -0.3e1/ , (16) 

where we have used the fact that the JT energetics for n=2, 5=1 is identical 
to that for n=l, 5=1/2 while for n=3, 5=3/2 there is no JT distortion. 
We have also used the t\ u orbital unimodal and bimodal splitting patterns of 
Part I to identify the filling {n\, n 2 , 72,3). In particular the fillings assumed are 
(0,0,2) for n=2, 5=0; (1,1,0) for n=2, 5=1; (0,1,2) for n=3, 5=1/2; (1,1,1) 
for n=3, 5=3/2. 

So long as Uw may be expected to vary slowly with the electron number 

n, then the two energy differences (|15j) and ( |i~6|) should be very similar. 

Moreover, the prevailing of a high- or of a low-spin state is decided by a 

fine balance between the Coulomb repulsion anisotropy (U\\ — U±) and the 

dynamical JT gain ^Eff s This suggests the possibility that if high-spin is 

more likely to prevail for C\q and C\q in the gas phase, where U is large, 

the balance might easily reverse in favor of low-spin when in matrix or in 

solution. Recent EPR data indicate that this is precisely the case. When 

frozen in a CH2CI2 glass, Cgo appears to be in a high-spin, 5=1 state [p0|] . 

Hence, in this case (U\\ — U±) is larger than 0.3 eV. However, optical and EPR 

data for Cqq in CH2CI2 and other matrices favor a low-spin state pi]. Now 
(3 s ) (3) 

Un — U]_ has therefore become smaller than 0.3 eV. We can conclude that, 
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even for a single embedded molecule, the balance between intra-ball Coulomb 
repulsion and dynamical JT energy gains is extremely critical. 

Recent photoemission and Auger data [33J have shown that the intra h u 
HOMO orbital Coulomb U is not as large as it was previously supposed. In 
particular, a decrease by only 0.23 eV from gas phase C 60 and the crystalline 
Ceo hole-hole Auger shifts, implies |U| < 1 eV in the latter. This upper bound 
is about a factor three smaller than those previously proposed ||35|| . In the 
light of this observation, it is not at all surprising to find that U\\ — U± is in 
the neighbourhood of 0.3 eV for Cqq in a matrix. 

Coming next to the pair binding energies of Ref 0, we find large negative 
dynamical JT-related U n 's for odd n. Even if we omit the A g contribution 
(which may be irrelevant for superconductivity, due to screening [|1|), we 
get U3=-0.2 eV. This negative value, will cancel at least a good fraction of 
the Coulomb positive intra-ball pair energy 

U£ oul = UP + (2f/{ 4) + 4f/j 4) ) - 2(t/{ 3) + 2C/| 3) ) « U\\ (17) 

This cancellation implies a severe decrease of the Coulomb pseudopotential 
H* relative to that calculated when the JT coupling is ignored fl3~6f . For a 
sufficiently strong solid-state screening of the electronic Un and U±, it may 
well be sufficient to reverse to a negative /i*, i.e. to an overall negative 
Hubbard U state. 

For n=2 and 4, dynamical JT stabilizes the average configuration of Cgg", 
since the pair energy is positive: U2 ~ 0.4 eV. This now acts to reinforce the 
bare Coulomb pair energy 

U^ oul = - 2U { 1 ] - 2Uf ] w 2U ± - U\\ 

jjCoul = ( 2[/ (5) + gt/ .(5)j + ^(3) + 2U (3)^ _ 2 ^ 2f/ (4) + AU (A)^ _ ^ _ ^ ^ 

where a filling (2,2,1) has been assumed for n=5. 

For even n, the JT coupling stabilizes a correlated insulating state of a 
lattice of evenly- charged C 60 molecules. In this type of insulator, fluctuations 
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about (rii) = n are suppressed, and a gap of order U n is opened in the 
electronic spectrum. 

This state has an even number of electrons per site, and is non-magnetic, 
very much like a regular band insulator. However, electron correlations re- 
sponsible for band narrowing and gap opening are vibronic in origin. We sug- 
gest that the (body-centered tetragonal f37|) structure of K4C60 and Pd^Cgo 
may be a realization of this state where electronic and vibronic interactions 



play an important role. So far, band calculations [23] and experiments [38 



had been in disagreement, the former suggesting a metal, and the latter 
finding a narrow-gap insulator. 



Very recent UPS data on K n C 60 JO have shown a decrease of the energy 



difference between the HOMO and the Fermi level (inside the t lu LUMO) 
when going from n=3 to n=4 and finally to n—6. This kind of non-rigid band 
behaviour is in itself not a surprising result. The surprise is that the decrease 
is very large from n=3 to n=4 (« 0.4 eV), and smaller from n=4 to n=6 
(« 0.2 eV). As pointed out by De Seta and Evangelisti, a positive Coulomb 
U would predict exactly the opposite. We observe that this behaviour is 
instead in agreement with our predicted pattern of effective U n of vibronic 
origin, which is therefore supported by these data. 

Additional experiments which may probe the electron-vibron interactions 



are short time resolved spectroscopy of excitons in neutral Ceo ||39|| . An 



exciton consists of an electron in the LUMO orbital and a hole in the H g 
HOMO levels, which interact with different strengths with the vibrons. The 
hole-vibron coupling inside the HOMO could be studied along similar lines 
to those presented above for the t\ u LUMO. 

As for superconductivity in solids with n=3, we expect the enhanced pair 
binding found here to be crucial for overcoming the on-site Coulomb repulsion 
and for enhancing T c over its value in, e.g., graphite intercalates. Broadening 
of the ti u electron level into a band of non-negligible width makes the present 
treatment insufficient for quantitative predictions. From the fundamental 
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point of view however, it is amusing to note that superconductivity can be 
enhanced by a decrease of lattice zero-point energy This adds to the usual 
BCS mechanism of reducing the electron kinetic energy by opening a gap. 
We hope to pursue this line of thought further in future work. 



6 Summary 

In conclusion, a full treatment of all the A g and H g modes has been given, 
and shown to yield analytical results with quantitative accuracy for the full 
dynamical JT problem of C^T. The ground state energetics has been studied, 
and unexpectedly large energy gains have been found, due to a decrease of 
zero-point energy. This implies large positive effective U n for n=2 and 4, and 
a large negative U3, which is very interesting in view of superconductivity in 
K3C60 and insulating behaviour in K 4 C6o- Detailed vibrational spectra for 
Cgo~ are presented, and proposed for spectroscopic investigation, particularly 
in gas phase. 

Related work is also being done by other groups pTT|. 
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TABLE I. Analytical expressions of energy shifts and excitation energies 
for the electron-vibron coupling of a single H g mode, for low-spin and high- 
spin states, to second order in the coupling constant g. 
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TABLE II. Vibronic excitation spectrum for the eight H g modes. Three 
different sets of coupling constants used in the perturbative expressions of 
Table I. The relations between the coupling strength g and the electron- 
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phonon coupling A/A r (ei?) (See Ref. are for H g modes g 2 = ^X/N(eF)/huj, 
and for A g modes g 2 = 3\/N(ep)/fiuj. 
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Table III 
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TABLE III. Analytical expressions for single mode pair energies (low-spin 
states) to second order in the corresponding coupling constants g^. 



26 



Table IV 





Exp. Energy 


coupling 


Ground state 




(cm^ 1 ) 


9k 


ener£ 


;y shift 


(meV) 




' 1 


[18| 


n—1 


n=2 


n=3 


E Ag,l 


493.0 


0.38 


-5 


-18 


-27 


E A g ,2 


1468.5 


0.39 


-14 


-54 


-81 


E Hg ,3 


270.0 


0.33 


-5 


-18 


-14 


E H g ,A 


430.5 


0.37 


-9 


-36 


-27 


E Hg ,5 


708.5 


0.20 


-5 


-18 


-14 


E Hg ,6 


772.5 


0.19 


-5 


-18 


-14 


E H g J 


1099.0 


0.16 


-5 


-18 


-14 


E Hg ,8 


1248.0 


0.25 


-12 


-48 


-36 


E Hg ,9 


1426.0 


0.37 


-30 


-120 


-90 


Eh 3 ,w 


1575.0 


0.37 


-33 


-132 


-99 


E a 






-18 


-72 


-108 


Eu 






-102 


-408 


-306 


E to t 






-120 


-480 


-414 


Un,A g 






-36 


18 


72 


Un,H g 






-204 


408 


-204 


U n jot 






-240 


426 


-132 



TABLE IV. Dynamical JT ground state energy shifts due to each mode, 
their total, and the pair energies U n . Results are accurate to second order in 
the coupling constants g^. 
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Figure 1: Exact (solid line) and perturbative energies (dashed line) for n—1, 
one coupled H g mode of frequency uo. The straight lines correspond to second 
order perturbation theory for the N v =0 and N v —1 multiplets. 



Figure 2: Exact (solid line) and perturbative energies (dashed line) for n—2, 
one coupled H g mode of frequency u. The perturbative lines correspond to 
the N v —0 and the lowest three multiplets of N v —1. 



Figure 3: Exact (solid line) and perturbative energies (dashed line) for n—2, 
one coupled H g mode of frequency u. The perturbative lines correspond 
to the N v =0 and the lowest three multiplets of N v —1. The second order 
coupling does not split the L=4 and L=2 levels in the N v —1 multiplet, while 
the exact theory finds they actually do separate. The next excited level (of 
species P) is also drawn to show its crossing of the initially lower D level 
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